Chiral Potts model and the discrete 
Sine-Gordon model at roots of unity 



Vladimir V. Bazhanov0 



Department of Theoretical Physics, 
Research School of Physical Sciences and Engineering, 
Australian National University, Canberra, ACT 0200, Australia. 



Abstract 

The discrete quantum Sine-Gordon model at roots of unity remarkably combines 
a classical integrable system with an integrable quantum spin system, whose param- 
eters obey classical equations of motion. We show that the fundamental i?-matrix 
of the model (which satisfies a difference property Yang-Baxter equation) naturally 
splits into a product of a singular "classical" part and a finite dimensional quantum 
part. The classical part of the i2-matrix itself satisfies the quantum Yang-Baxter 
equation, and therefore can be factored out producing, however, a certain "twist" of 
the quantum part. We show that the resulting equation exactly coincides with the 
star-triangle relation of the A^-state chiral Potts model. The associated spin model 
on the whole lattice is, in fact, more general than the chiral Potts and reduces to 
the latter only for the simplest (constant) classical background. In a general case 
the model is inhomogeneous: its Boltzmann weights are determined by non-trivial 
background solutions of the equations of motion of the classical discrete sine-Gordon 
model. 
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1 Introduction 



The fundamental works of Professor Akihiro Tsuchiya [1, 2J made an outstanding contri- 
bution to the theory of integrable quantum systems. Here we present some new devel- 
opments in this field. The discovery of the chiral Potts model originated in [3l SI [5] and 
finalized in [6] brought a lot of new interesting problems. First the Boltzmann weights 
of the model require high genus algebraic functions for their uniformization. The second 
(related) problem is that there is no difference property in the model. The uniformization 
problem has been solved in [7] using classical results on the rotation of a rigid body [8]. 
However, even with this uniformization we still do not know how to apply various methods 
which worked perfectly for all other two-dimensional solvable models, basically because 
there is no a difference property. In particular, it took 16 years until the intriguingly 
simple conjecture of ^ for the local order parameter 

< 0.°^- >= (A;Y^ ^^ = ^-^^^ j = l,...,iV-l. (1.1) 

has finally been proven [10] by a generalization of the method of [11]. Here N, N >2, is 
the number of local spin states and is a modulus of the the algebraic curve appearing 
in the model, it also plays the role of the temperature like variable: the model is critical 
for = 0, it reduces than to the Fateev-Zamolodchikov model |12]. For N = 2 the 
chiral Potts model coincides with the ordinary Ising model and formula (11. ip reduces to 
the famous Onsager-Yang result |13l Hlj for the spontaneous magnetization of the Ising 
model. 

In view of the above discussion the chiral Potts model may seem standing very far 
away from the "conventional" solvable model with the difference property. However, this 
is not so. It was shown [15] that the chiral Potts model is remarkably connected to the 
six-vertex model: its R-matrix satisfies the Yang-Baxter equations with two cyclic L- 
operators related to the R-matrix of the six-vertex model. Further, it turned out that the 
last connection can be extended by replacing the six-vertex model (with two states per 
edge) by an the n-state model [12] associated with the Uq{sl{n)) algebra. This resulted 
in the "s/(n) -generalized chiral Potts model" [1^, which is a two-dimensional model with 
spins that each take A^"~^ values. 

Note, also that the chiral Potts model has very deep relations with the theory of 
cyclic representations of quantum groups at roots of unity [IHl [Hj- In particular, the 
R-matrix of the model can be interpreted [171 [IH] as an intertwiner for two minimal cyclic 

representations of the Uq{sl{nj) algebra. 

Next, in [20] it was shown that the s/(n) -generalized chiral Potts model can be in- 
terpreted as a model on a three-dimensional simple cubic lattice consisting of n square- 
lattice layers. At each site there is an A^-valued spin. And when N = 2 this three- 
dimensional model reduces (to within a minor modification of the boundary conditions) 
to the Zamolodchikov model f2T[ [22] . For > 2 it gives a new solvable interaction- 
round-a-cube model in three dimensions [20]. The Boltzmann weights of this model 
satisfy [231 1211 the tetrahedron relations which ensure its integrability. However, it turns 
that the integrability of the model can be proved bypassing the tetrahedron relations. It 
is based on a much simpler "restricted star-triangle relation" [25J, which is a particular 
case of the star-triangle relation [SI ES] for the original chiral Potts model. 

In this paper we discuss yet another remarkable connection of the chiral Potts model. 
This time we relate it with the quantum discrete sine-Gordon (QDSG) model [2711281 129]. 
In [30j it was shown that for rational values of the coupling constant the QDSG-model can 
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be viewed as some special two-dimensional lattice model of statistical mechanics where 
(integer) spins take N > 2 values and where the Boltzmann weights are determined 
by solutions of the equation of motion of the classical discrete sine-Gordon model (the 
situation here is very similar to quantum field theory in a classical "background" field). 
Here we show that in the simplest case of the constant background the above lattice model 
exactly coincides with the chiral Potts model. 

The organization of the paper is as follows. In Section 2 we review the formulation and 
basic properties of the classical and quantum discrete sine-Gordon models. Associated 
solutions of the Yang-Baxter equation are considered in Section 3. In Section 4 we review 
the basic definitions of the chiral Potts model and Section 5 we relate this model to the 
QDSG model. 



2 The quantum discrete Sine-Gordon model. 
2.1 Formulation of the model. 

The Sine-Gordon equation for a scalar function 0(x, t) in 1 + 1 dimensions 

- + (9^0 = sin (2.1) 

is the most famous example of a system integrable on both the classical and quantum 
levels. Despite of an extensive literature devoted to this model (see [31],[32] for the 
references), it still continues to reveal its new features. In this paper we shall consider 
integrable generalizations [271 [281 [29] of this model (classical and quantum) for the case 
when both the space and time variables are discrete. 

Following [29] let us give a brief description of the discrete quantum Sine-Gordon 
model. The classical model can then be obtained in an appropriate limit. Let q be a 
complex number, and Aii^) be an algebra of power serie^ in variables Wn, n = 0, . . . ,2L — 
1, L > 2, obeying the following relations 

WnWn^l = q'^Wn-lWn, Vn, 

[2.2) 

WnWm = WmWn, \m - n\ > 2, 

with the periodicity condition w„+2L = Wn- The variables w constitute the set of the 
dynamical variables of the model at any fixed value of time. The evolution on one step 
of the discrete time acts as an automorphism r of the algebra Aiiq) 

r : A(q) — . A(q), (2.3) 

such that 

U!2n =r{w2n) = /(qW2n--l) ^2™ /(qW2n+l)"\ 

(2.4j 

U)2n-1 =r{w2n^l) = /(qW2n-2) W2„_l /(qt&2n)" • 

Here and below we use the notation f{x)^^ = l//(x). The function / has the form 

1 -I- Ar 

/(x) = /(K^x), f{X,x) = ^^ (2.5) 



^ We shall assume these series to be semi-infinite, i.e, the powers of w to be bounded from below, but 
not necessarily non-negative. The multiplication of such series is well defined. 
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where is a fixed parameter of tfie model. One can easily show that the variables Wn 
satisfy the same relations (12. 2p as Wn- To get a geometric picture of the equations (12.41) 
consider the square lattice drawn diagonally as in Fig. [1], with L > 2 sites per row and 
impose periodic boundary conditions in the horizontal (spatial) direction. The time axis 
is assumed to be directed upwards. Now assign the variables w„, n = 0, . . . , 2L — 1, to 
the sites belonging to some horizontal "saw" S as shown in Fig 1. In the same way assign 
the variables tin to the saw S shifted by one time unit from the initial saw S. 



W2n W2n+2 



W2n-1 W2n+1 

W2n W2n+2 



S 



Figure 1 - Arrangement of the operators Wn and Wn on the space-time lattice. 

With such arrangement both evolution equations (12.41) connect four w's around a single 
square face of the lattice and can be written in a universal form 



(2.6) 



where the variables w are labeled as in Fig. [2J 



^UJD 



Figure 2 - Labeling of the four w^s around an elementary face. The indices stand for "up" , 
"down", "left" and "right". 

The algebra (12. 2p has two Casimir elements 

L-1 L-1 

Cl = n ^2n+l , C2 = X[w2n. (2.7) 

?i=0 n=0 

which commute with all w's and are preserved by the evolution automorphism r. Here 
we restrict ourselves to the case 

Ci = C2 = l. (2.8) 
This leaves 2L — 2 independent generators, say, wq, wi, . . . , W2l-3- 
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The quasiclassical limit of the model corresponds to the case q —>■ 1. The algebra (12.21) 
is then replaced by the Poisson algebra, where 

{Wn,Wn-l} = 2WnWn-l, Vn, (2.9) 

but all other brackets {wn,Wm} vanish. The equation of motion (12. 6p becomes 

Wu = WDfiwL)/fiwR). (2.10) 

One can show that in the appropriate continuous limit the equation (I2.10p reduces to 
(EUD when 0. For the details see 

We will say that the evolution r is Hamiltonian if there exists an invertible operator 
U such that 

T{a)=V-^aV. (2.11) 
Let r(A,x), A G C, be a solution of the following q-difference equation 

r{X,q^x) = f{X,qx)r{X,x), (2.12) 

where the function / is given by (12.51) . Define the operator U as follows 

L-l L-1 

V = l[r{K\w2k)l[r{K\w2k+i). (2.13) 

fc=0 fc=0 

Applying (12.111) for a = Wn and using relation (12.121) one can easily see that (12.111) reduces 
then to the evolution equations (12.41) . The calculations are rather trivial due to the locality 
of the commutation relations among w's. 

To proceed further we need to specify an actual realization of the Hilbert space of 
the model, supporting the commutation relations (12. 2p . where the equation (I2.12p is well 
defined and has a unique solution. At the moment there is no classification of representa- 
tions of this algebra suitable for this purpose. Only two cases are known so far. The first 
case is connected with non-compact modular representations of (12. 2p and the Faddeev- 
Volkov model [271 129] (see also [33l [Ml [35l [36]). Recently [37i i38i| this model was solved 
exactly and shown to be related to the discrete analog of the Riemann mapping theorem 
in two dimensions. For further details we refer the reader to the already mentioned papers 
[371 [38] and to [39]. 

In the present paper we consider the second known case when the equation (I2.12p is 
well defined. It is when the parameter q is a root of unity. Originally this case was studied 
in m- 

2.2 The model at roots of unity. 

Here we are continuing the study of the discrete quantum Sine-Gordon model in the limit 
when q^ approaches a root of unity. As it was shown in [30] the model in this case can 
be viewed as a (finite dimensional) integrable quantum system on top of an integrable 
classical one. 

Let (— qo) be a primitive A^-th root of —1, then c\q will always be a primitive A^-th 
root of 1, 

(-qo)^ = -l, qf = l, Ar>L (2.14) 
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When q = qo the algebra ^^(qo) has a (commutative) center Z{ALi^o)) generated 
by the elements . It is easy to obtain the evolution equations for these commuting 
variables. Taking the A^-th power of both sides of (12. 6p for q = qo one obtains 

w^ = w^f{^'\w^)/f{n'\w^), (2.15) 

by using (12.141) and the simple identity 

N~l 

/(^^^x^) = ^/('^^qr^), (2.16) 

j=0 

where f{\,x) is defined in (12.51) . Apart from the trivial replacement of by the 
A'"-th powers of w obey the same evolution equation as the one of the classical discrete 
Sine-Gordon model (I2.10p . In particular, when = 1, q = 1, formula (12.151) exactly 
reduces to (12.101) as it, of course, should. 

The quantum evolution of the model takes place in a finite dimensional factor of the 
algebra ^L(qo) over its center Z{Al{c{o))- In fact, let a = {oq, . . . ,a2L-3} be a set of 
nonzero complex numbers and la be an ideal generated by {w^ — an), n = 0, . . . ,2L — 
3. Then the factor AL{c\o)/Ia for any a is isomorphic to a finite dimensional algebra 
generated by elements Xi, Zi, i = 1, . . . , L — 1 obeying the following relations 



XiZi — c\q ZiXi, — Z^ 

XiZj = ZjXi, i ^ j. 

The explicit isomorphism is achieved by the formulae 



(2.17) 



W2n = e^"+^Xn+i, n = 0, . . . , L - 2, 

W2n+i = e«"+^-«"+^Z-^iZ„+2, n = 0, . . . , L - 3, (2.18) 

W2L-3 = e'''^-'Zj^\ , 

where for later convenience we have parametrized the a's through a set of new variables 
Pi) Qi-i i — 1,...,-Z^ — 1 as 

a2n = e^"+i, n = 0,...,L-2, 



n = 0,...,L-3, (2.19) 

a2L-3 = e^' 



The remaining two dependent generators W2L-2 and W2L-1 can be expressed from (I2.18p . 
(1X71) and (12:81) . 

One can introduce a Poisson algebra structure on Z{Al{(\o)) with its Poisson action 
by derivations on the whole ^^(qo)- In fact, this structure naturally emerges from (12.21) in 
the limit q — >^ qo. First, recall that the Poisson algebra P is by definition a commutative 
algebra with the structure of a Lie algebra, defined by a bracket { , } and the rule 

{a, he} = {a, b}c + b{a, c}, (2.20) 

for any a,b,c G P, which determines a distributive action of the bracket on the (com- 
mutative) product. Next, we can identify the algebras Al{^) for all q as vector spaces 
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of formal power series in generators by choosing some definite way of their ordering. For 
convenience set 

q = e^o (2.21) 

where h is a formal variable. Denoting by a ■ 6 the product of two elements in .4.L(q) we 
can define the bracket 

, ^ , a ■ b — b ■ a 
|a, b] = lim 2.22 

h^O h 

when at least one of the elements a and b belongs to Z{Al{c{o))- Using (12.21) it is easy to 
compute the Poisson brackets (12.221) between the generators of ^(^^(qo)) and their 
action on 



(2.23) 



where 



1, 171 = 11 (mod 2L), 

^ ' (2.24) 

0, otherwise. 

The Poisson action of Z{Al{<^o)) on the whole Al{c\o) is extended from (12.231) by the same 
rule (12.201) where now a G ^(^^(qo)) and 6, c G ^^(qo) (of course, the multiplication there 
should now be assumed non-commutative). 

For later use introduce the following notation 

exp{a)ob = ^—^{a---{a,b}---} (2.25) 

for the Poisson action of exp(a) on b, where a G Z{Al{<^o)) and b G v4L(qo). It follows 
from (12:221) that 

exp (j^) - b- exp (^-^) = exp(a) ob + 0{h) (2.26) 

If we use the parametrization (I2.18P then (I2.23P reduces to the only non-trivial brackets 

{Q^,P,} = 6,, (2.27) 

The standard quantization of the Poisson structure described above gives a realization 
of the algebra ^^.(q), q = e'^qo, within a wider algebra of formal power series 

X{P, Q, X, Z, h) = (series in Pi,Q„X„Zi, h) (2.28) 

in non- negative powers of Pi, Qi, Xi, Zi, i = 0, . . . , L — 1 and h, where Xj, Zi are defined 
in (I2.17P and Pi, Qi satisfy the relations 

[Q„P,] = 2M.„ [P,,h] = [Q„h] = 0, ^^^^^ 
[P„ X,] = [P„ Z,] = [Q„ X,] = [Q„ Z^] = 0. 

The generators w„ are given by (12.180 where the elements Pi,Qi now satisfy (I2.29p . This 
realization of Al{c\) is most suitable for the study of the limit q — > qo, qo^ = 1. 

Note also that when L = 2 the algebra Al{c\) has only two independent generators wo 
and Wi which we denote as v and u respectively. From (12. 2p we have 

uv = q^vu (2.30) 
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while the reahzation fl2.18p contains only one pair of the P, Q and one pair of the X, Z 
operators 

u = SZ-^ V = e^X (2.31) 

where 

XZ = ql ZX, X^ = Z^ = 1, (2.32) 

[Q, P] = 2h [P, h] = [Q, h] = 0, (2.33) 
where h and q are related by (12.211) . 

2.3 The integrability of the model. 

In this section we briefly address the matters related to the integrability of the discrete 
quantum Sine-Gordon model which so far have not been discussed. The details of the 
calculations can be found in [3Qj- 

Let Wq be the Weyl algebra generated by the invertible elements U and V, 

Wq : UV = qVU , 

where q G C. Further, let -R(A) be the i?-matrix of the 6-vertex model acting in (8> 

R{X) = (Aq - A^^q-^) (en ® en + 622 ® 622) 

+ (A-A-i)(en®e22 + e22®en) (2-34) 
+ (q - q"^)(ei2 ® e2i + e2i ® ei2) , 

where Cij is the matrix unit and L{X) be an L-operator acting in Wq 

^(^' = (aK-' (2,35) 
These operators satisfy two Yang-Baxter equations 

i?12(A)i?i3(AAi)i?23(/i) = i?23(/x)/2l3(A/i)i?12(A) (2.36) 

and 

Ri2i\)Lii\fi)L2ifx) = L2(^)Li(A/x)i?i2(A) (2.37) 
which are written in the standard notations. As usual, define the transfer matrix 

t{X) = trace(Lo(AK)Li(A//«)L2(AK) ■ ■ • L2l-i(A/k)), (2.38) 

for a periodic chain of the length 2L with alternating rapidity variables Ak, A/k, Xk, etc., 
where k is a constant. This transfer matrix acts in a direct product \N^^^ of 2L copies of 
the algebra Wq. 

As a consequence of (12.371) the transfer matrices (12.51) with different A form a commut- 
ing family. From the definitions (I2.35P and (I2.38P it is trivial to see that t{X) is a polyno- 
mial of the degree 2L in A^. What is less trivial, the operators Ui, Vi, i = 0, . . . ,2L — 1, 
associated with different sites of the chain, enter coefficients of this polynomial only trough 
2L elements of the form (remind that we work with periodic boundary conditions) 

Wn = UnV-'Un+lVn+l, (2.39) 
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which obey the commutation relations of the algebra (12.21) ! The proof of this statement can 
be found in [30]. The next remarkable fact is that the transfer matrix (12.381) commutes 
with the evolution operator (12.111) and, therefore, it generates commuting integrals of 
motion of the discrete quantum Sine-Gordon model. 

To see this let us search for another i?-matrix acting in Wq ® Wq which intertwines 
two L operators in the "quantum" space. In fact, this i?-matrix has already appeared in 
Sect. 1. Assume that there exists a well defined element r(A,M) G Wq Wq, where 

u = UV-^ ® UV, (2.40) 

and r(A,x) is the solution of the q-difference equation (12.121) . 

r(A,q2x) = \±^r(A,x), (2.41) 
A + qx 

Then, it is easy to show that 

r(A, M)(L(A/i) ® Lifi)) = (L(/i) ® L(A/i))r(A, u) (2.42) 

With an account of (12.111) . (12.391) and (12.401) this fact immediately implies the required 
commutativity of the evolution operator of (12. lip and the transfer matrix (I2.38P 

[t(A),U]=0 (2.43) 

At this point we would like to note that so far we used only the defining relation (12.411) 
for the i?-matrix r{X,x) rather than its explicit form. This concerns the Yang-Baxter 
equation (12.421) and the equivalence of the two forms (12. 4p and (12. lip of the equations of 
motion in the model. An explicit form of r{X,x) for the case when q^ is approaching a 
root of unity will be given in the next Section. 



3 The Yang-Baxter equation for r(A, x). 
3.1 General analysis. 

The reader may have noticed that the sequence of Yang-Baxter equations (12.360 . (I2.37p . 
(12.420 in the previous Section which starts from the six-vertex model is very much resem- 
bling the one that appeared in the analysis of the chiral Potts model in ref. [15]. In fact, 
our considerations here are quite similar but more general than in [15] and the relations 
with the Chiral Potts model will soon be very explicit. 

The final point in the descending sequence of the Yang-Baxter equations mentioned 
above is the equation for i?-matrix r(A, x) itself. The form of (12.420 requires that it should 
be written as the "braid relation" 

r(A, u)r{Xfi, v)r{fi, u) = r{fi, v)r{Xfi, u)r{X, v), (3.1) 

where v and u denote respectively any two successive generators Wn and w^+i from (12.400 
which form a Weyl pair 

uv = q^vu. (3.2) 

Here A, G C are the (multiplicative) rapidity variables. It might be tempting to assume 
that the two relations (I2.4ip and (13. 2p are sufficient to imply the Yang-Baxter equation 
(13. ip . As we shall see below this statement is not in general true. However, it would be 
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quite interesting to understand how far one could advance in proving (13.11) using those 
two relation only. 

Combining (12.411) and (13. 2p one obtains 

r(A, u)v = V ^ — r(A,M), r(X,v)u = u- — ^—r(X,v). (3.3) 
X + qu q + Af 

Denote the left and right hand sides of the Yang-Baxter equation (13.11) as $l and $r 
respectively. Then using the relations (I3.2p . (13.31) only one can show that for both $ = $i 
or $ = 

= ug{u,v,q)^, = v(^{g{v,u,q~^))~^, (3.4) 

where 

g{u, V, q) = (qA^/i + q^Xfiu + v + qXvu){qX + q^u + Xfiv + X"^ fivu)^^ . (3.5) 

The calculations are elementary but a bit tedious. We present them in Appendix A. 
Fortunately, we have to derive only one of the equations in (13. 4p since they follow from 
one another under the action of the automorphism 

r(A,x) — >{r{X,x))~^, u< — > v, q — ^q~^, (3.6) 

which leaves (I3.2p . (l3.3p unchanged. There are also additional "exchange" relations which 
follow from (13.40 and the fact that r{X~^,x)~^ satisfies the same q-difference equation 
fICTD as r(A,x) 

= u^{g{u,v,q))~^, ^ v = vg{v,u, q~^)^. (3.7) 

As above, these relations are valid for both $ = $i or $ = The function 'g{u,v,q) 
is given by the RHS of (13.50 with A and /i replaced by and A~^ respectively. The 
relations (13. 4p and (13.70 immediately imply that 

[^,u] = [^,v] = 0. (3.8) 

where 

^ = ^L^R (3.9) 

Apparently, this is the most general consequence for the Yang-Baxter equation (13.10 
which one can obtain from (13. 2p . (13. 30 without further requirements on the properties of 
the r(A,x). Therefore, if the representation of (12.20 were such that the commutativity 
with the elements u and v in (13. 8p implies that the quantity \1/ is proportional to the 
unit operator, then Yang-Baxter equation (l3.ip would be satisfied and, at the same time, 
equation (12.120 would have a unique solution. 

Let q = e'^qo, where h is a formal variable and qo is defined in (12.140 . The most 
appropriate way to study (I2.12p when q^ approaches a root of unity, /i ^ 0, is to consider 
r(A, x) and, hence, the left and right hand sides of equation (13. ip as asymptotic expansions 
in h. Define the Euler dilogarithm function as 

U,{x) = - T-^^^dt. (3.10) 
Jo 't 

The Yang-Baxter equation (13. ip involves only one Weyl pair (13.20 which we parametrize 
as in (I2.3ip . Then it is convenient to rewrite (I2.4ip as 

r(A,e^+'^qo^) = ^^r(A, e^-'^z/qo). (3.11) 
A + a; 
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Here we set x = e^z, where z is a formal variable such that z'^ = 1. Next, impose the 
following normalization condition on r(A,x) 

r(l,x) = l. (3.12) 

Then equation fl2.4ip uniquely determine the asymptotics expansion of r(A, x) in the form 

r(A, e^z) = exp {c_i(T, z)h-^ + co(T, z) + Ci{T, z)h + ...}, (3.13) 

where the coefficients Ck{T, z) are series in T and z. In principle, all these coefficients can 
be calculated explicitly from (13.111) . In particular, retaining only the first two terms in 
(13.131) one obtains 

r(A,x) = exp ^ ^ r(A, x)(l + 0(/i)) (3.14) 



where x = e^z. 



and 



H{a, h) = {Li2{-ab) + Li^i-a/b) + log^ b + vrVe} (3.15) 

/ \N \ {N-1)/2N N-1 / x2j+l \ 

The branches of the function Li2{x) are assumed to be chosen such that H{l,x) = 0, so 
that the normalization condition (I3.12p is satisfied. The expansion (I3.14p can be easily 
deduced from formula (3.11) of ref. [40j . 



3.2 The case q ^ 1. 

First, consider the case q — > 1. This corresponds to = 1, qo = 1, q = e'^ in our previous 
notations. Eq. (l2.3ip now reads 

u = e^, v = e^. (3.17) 
It follows from (I3.13P that both elements r(A, u) and r(A, v) have the form of exponentials 

exp{X/h), (3.18) 

where X is a. series 

oo 

^=5^ CMmh'Q'P"', (3.19) 
k,e,m=0 

in non-negative powers of P, Q and h. Of course, for r(A, u) and r(A, v) the above series do 
not contain mixed terms (they involve only pure powers of P or Q). However, below we will 
consider more general series fl3.19p . We will always assume them to be "normally ordered" 
such that all P's appear to the right from Q's. The product of such two exponentials is 
again an exponential of the same form. Indeed, from the Campbell-Hausdorff formula we 

e'^/^y/'^ = e^/^ (3.20) 
^ = + ^ + 2^^-^' ^] + ^f'^' f'^' + '^l]} + • • ■ (3.21) 
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where for each term of the last series the number of commutators coincides with the power 
of the parameter h in the denominator. If both X and y are of the form f l3.19p . then 
due to f l2.33p each commutator produces one extra power of h in the numerator so that 
all negative powers of h in (13.211) exactly cancel out. Therefore Z is a. series of the form 
(13.191) as well. It should be noted that each coefficient therein will now be an infinite 
numerical series in the coefficients of X and 3^. For the moment let us assume that all 
these numerical series converge and, hence, Z is well defined. Repeatedly applying (13.201) 
one concludes that the ratio \1' of the left and right hand sides of the Yang-Baxter equation 
in (13. 9p is also an exponential of the form (I3.18p . Then it is easy to show that the most 
general form of ^ satisfying (13. 8p . (13.170 is 



= const exp {i7r(mP + , (3.22) 



where m and n are arbitrary integer constants. Finally setting the spectral parameters 
A = /i = 1 in (13. ip and using (I3.12p one concludes that = 1. Since we worked with 
exponentials of asymptotic series in h this implies that the Yang-Baxter equation (13.10 is 
satisfied as an equality of such exponentials in any given order in h. 

It should be noted that the series in h in the exponents of (13.130 and (13.200 can be 
consistently truncated at any level. In particular, setting N = 1 and taking into account 
only the two leading terms in (I3.13P one obtains the following Yang-Baxter equation 

r(A, u)r{Xfi, v)r{fj,, u) = r{fj,, v)r{Xfi, u)r{X, v){l + 0(h)), (3.23) 

where 

r(A, x) = e^(^'")/'^(l + Oih)). (3.24) 
Writing (I3.23P in full and substituting (13.170 one obtains 



exp{H{X,e^)/h) exp{H{Xfi,e^)/h) exp{H{fi,e^)/h) 

= exp{H{ij,eP)/h) exp{H{Xi2,e'^)/h) exp(if(A, e^)//i)(l + 0(/i)) 



(3.25) 



Introducing a special notation for the leading term of the Campbell-Hausdorff composition 
in fl3:20|) 

exp (^j^ exp (^^^ = 0(1) exp , h ^ 0, (3.26) 

we can rewrite the equality of the singular-in-/;, exponents in (I3.23P as 

H{X, u) * H{Xii, v) * H{fi, v) = H{fi, v) * H{Xfi, u) * H{X, u). (3.27) 
Below we show that this formula is equivalent to some identity for the dilogarithm function 

(EinD. 

To complete the proof of the Yang-Baxter equation we have to ensure that the above- 
mentioned numerical series involved in the product of the i?-matrices converge. Appar- 
ently this can be done "order by order" in perturbations in h using a special structure 
of the Campbell-Hausdorf series (I3.2ip and the fact that the numerical coefficients there 
rapidly decrease (see, e.g., [41J, chapter 2, §6.4). We will not analyze this problem in its 
general setting here, but restrict our considerations to the truncated equation (I3.23p . In 
this case we are able not only to show that products of the /^-matrices in (I3.23P exist 
but also explicitly calculate normal symbols of the two sides of (I3.23P up to the order of 
0{h^). The details of the calculations are presented below. 
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3.3 The 12-term dilogarithm identity 

Let A be an operator given as a power series in P and Q satisfying f l2.33p 

oo 

A=J2 ""^iQ'^' (3.28) 

k,l=0 

ordered such that all P's are placed on the right from Q^s. We will call the series ordered 
in this way as normal form of the operator A. As usual the normal symbol : A{Q, P) : of 
the operator A is defined by the same series fl3.28p as its normal form but with P and Q 
replaced by commuting variables P and Q respectively. 

To calculate the normal symbol of the product of two different operators one has to 
use a simple identity for the normal symbol of an elementary (unordered) monomial 

Am = P'Q' (3.29) 

: AM,P) := (-2h)''+' exp{PQ/2h)^-^expi-PQ/2h) (3.30) 

which can be easily proved from (12.331) by induction. Following [12| we can write the RHS 
of (I3.30p in an integral form 

= / exp((Q - Q')iP - P')/2h)P'^Q'' (3.31) 

J 2TTh 

assuming h to be negative real and considering P' and Q' as complex conjugated variables 
P' = [Q'Y G C. The integral is assumed to be over the whole complex plane. Then the 
normal symbol of the product of two operators A and B can be written in a compact way 

f dO'dP' 

■.ABiQ,P):= j ^±-^expiiQ - Q')iP - P')/2h) : AiQ,P') B{Q',P) : (3.32) 
Applying this for the left hand side of (I3.23P one obtains 

: $i(Q, P) := r(A, e^^) J ^iQ^Q')iP~ni^h ,(A;,, e^') r(/., e'^') (3.33) 

This integral can be explicitly calculated to within terms of the order of 0{hP) using the 
saddle point approximation. A similar treatment of the quantum five-term identity can 
be found in [13]. Substituting fl3.24p into (13.330 one obtains 

: $l(Q,P) := {detyi{x\y'))-"^eMEL/2h){l+0{h)) , ^^^^^ 
El = H{\ x) + H{\^, y') + i7(/i, x') + 2 log(x/x') \og{yly') 

where 

a; = e^, y = e^, x' = e<3', = e^', (3.35) 

and the coordinates P', Q' of the saddle point are determined by the equations 

^^A^VMl y^y:^l±t^, (3.36) 

Xjji + y' 11 + x' 



^ In this subsection we will denote the operators P and Q from (|2.33p as P and Q and use the letters 
P and Q for commuting arguments of the normal symbols. We hope this will not cause any confusion. 
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The two by two matrix M coincides (up two a factor of 2h) with the matrix of the 
quadratic form in the exponent of the integrand near the saddle point 



M{x',y' 



( x' (/x^ - 1) ^ \ 

(x' + /i) (1 + 

, y' (A V - 1) 

\ (l + A/i?/0(?/' + A^)/ 



(3.37) 



The equations 03.361) determine two distinct saddle points. Instead of solving these equa- 
tions with respect to x' and y' , it is more convenient to regard x' and y' as independent 
variable and consider (13.361) as a definition of x and y. In the same way one can calculate 
the normal symbol of the RHS of (I3.23p . Then for any x' and y' the corresponding saddle 
point for : ^r{Q, P) : is given by 

yi, ^ x\\nx'y' + x' + \y' + ji) ^„ ^ y'jXfix'y' + x' + Xy' + AV) 
fix'y' + Xx' + + A/i ' X'^fix'y' + Ax' + ?/' + A/i 

where 

a:" = e^", / = e«", (3.39) 

and P" and Q" are integration variables in the corresponding integral for : ^ji{Q,P) :. 
The expression for : $/j(Q,P) : is obtained from (13.341) merely by replacing x' and y' 
therein with x" and y" respectively. Equating the singular-in-/i exponents of the normal 
symbols we thus obtain the following identity. 

The twelve-term dilogarithm identity. Let A, /i, x' and y' be arbitrary complex numbers 
and X, y, x", y" are given by (I3.36p . (13.381) . Then 

H{X, x) + H{X^i, y') + iJ(/i, x') + 2 \og{x/x') \og{y/y') 

(3.40) 

= i/(/i, x") + H{X^i, y") + H{X, y) + 2 log(x/y") \og{y/x"), 

where the function H is defined by (I3.15p . The branches of the dilogarithms and loga- 
rithms are to be chosen as analytic continuation from the case when A, /x, x' and y' are 
all real and positive. 

Since we have shown (by explicit calculations) that products of the i?-matrices in the 
Yang-Baxter equation (I3.23P do exist and are well defined (to within the required order 
in h) the identity (I3.40p is a corollary of our proof of the Yang-Baxter equation in the 
previous subsection. However, is very easy to check directly the dilogarithm identity of 
this type once it is written. First, computing partial derivatives of the first order in the 
variables A, /x, x' or y' one can check that the difference between the LHS and RHS of 
(I3.40p considered as a function of these four independent variables is a constant. To check 
that this constant is actually equal to zero one has to set x' = y' = 1. The four variable 
dilogarithm identity (13.400 seems to be new, but, of course, can be obtained by a repeated 
application of the classical five term identity for the dilogarithm function (I3.10p (for the 
most comprehensive survey on this subject see [Hj). 

Note that ^m^ . flX^ imply 

..?Al±i^. ,^£!(i±MO. (3.41) 

/i -|- x" A/x -|- y" 
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Regarding A and /i as constants and x", y" as independent variables we can write the 
RHS of f l3.40p as some function J^o{x",y"). Then the identity (13.401) can be interpreted as 
an invariance of this function under the (bi)rational substitution f l3.38p of two variables 

W,y") = W,y') (3.42) 

Note that the substitution (13.381) is an involution, i.e, it reduces to identity substitution 
on the second iteration. 

Similarly, the order of 0{h^) in (13.231) (which is equivalent to to the equality of the 
determinants in (I3.34p and in the corresponding expression for : $ij(Q,P) :) gives a 
polynomial invariant 

jr^(x", y") = (1 + i2x"){l + Xfiy") (3.43) 

of the substitution (13.380 

J^,{x",y")=J^,ix',y'). (3.44) 
Of course, the last formula can be readily checked by direct calculations. 

3.4 The case of a generic root of unity 

Now let qo be a primitive A^-th root of —1 as defined in (12.140 and 

q = e^o (3.45) 

The generators u and v are now given by the general formulae (I2.3m2.33|) and the leading 
asymptotics of the i?-matrix by (13.141) . Our previous proof of the truncated Yang-Baxter 
equation (13.231) can be easily modified for this case. One has to generalize the series in 
(13.191) to series in P, Q, X, Z and h (of course the powers of X and Z will not exceed 

— 1). The subsequent arguments are very similar to those of the case q — > 1 and we 
will not present them here. 

Let us rather consider consequences of the Yang-Baxter equation (I3.23P in this case. 
The i?-matrix (13.140 factorizes into a "classical" and a "quantum" part. The former is the 
singular-in-/i exponential belonging to the center of the algebra Z{A2{c\o)) while the latter 
belongs to the finite dimensional factor of this algebra over its center, which is generated 
only by the elements X and Z obeying (12.320 . It turns out that that the Yang-Baxter 
equation (I3.23P remarkably splits into two separate (quantum) Yang-Baxter equations for 
the classical and the quantum parts of the i?-matrix. The trick is similar to that used in 
[IQ] for the asymptotics of the quantum five-term identity of ref. [13]. The equation for 
the classical parts is equivalent to (13.251) while the one for the quantum parts appears to 
be equivalent (as we shall see in the following sections) to the star-triangle relation of the 
A^-state chiral Potts model. 

Substitute (I3.14p into (I3.23P and move all singular-in-/i exponentials to the right. To 
do this one has to use formula (12.261) . since these exponentials in general do not commute 
with the elements r. For instance, for the LHS of (13.230 one obtains 

^^(A,e«)//i ^x,e^Z-') e^(V,e^)//^ f(A/x,e^X) e^(/^' e'^)/^f (/x,e«Z-i) = 
r(A, e«Z-i) (e^(^' e"^) o r(A/i, e^X)) (e^( V, e^) o e^(/^' ^'^) o r(;x, e^Z-^)) x 
^,#(A, eQ)/h^H{\^^, e^) ^^^(/x, e«)//.(i + o{h)), 

(3.46) 
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where 

H{a,b) = N-^H{a^,b^) (3.47) 

with the function H given by fl3.15p . Perform similar transformations for the RHS of 
(13.231) . After that all singular exponentials in the LHS cancel with those in the RHS side 
due to fl3.25p (more precisely, one has to use (13.251) with the function H replaced by a 
function H, defined in (13.471) ). The remaining equation is that for the quantities of the 
order of 0{h^) so that one can now set h = 0. The operators P and Q then become 
commuting variables. More precisely, they become elements of the Poisson algebra with 
the bracket 

{Q,P} = 2h (3.48) 

as it was explained in Section 2.2. As a result we obtain the following "twisted" Yang- 
Baxter equation for the elements (13.16^ 

r(A, e«Z-i)r(A/i, e^'A:)r(/i, e'^" Z'') = r(/i, e^X)r(A/i, e«'Z-i)r(A, e^"X) (3.49) 

Here 



P' = exp (^H{X, e^)^ o P, 
g' = exp (^^(/i,e^)) oQ, 
P" = exp (^H{X, e^)^ o exp (^^(A/x, e^)^ o P 
Q" = exp (^H{fi, e^)^ o exp (^^(A/x, e^)) o Q . 
where we have used the notation (I2.25p . Using (I3.48P one can show that 

]\fp' i-r/v e NP" i-|-A/ic 

^ X^ + e^Q ^ ' ^ A^^ + e^Q' ^ ' 

NQ' ^ ^ NQ NQ" ^ ^ ^ NQ 

1 + fi^e^P e , e ^ ^ X^fi^e^P' 



(3.50) 



(3.51) 



The Yang-Baxter equation (13.490 contains four (complex) parameters: two rapidities 
A and and two arbitrary parameters P and Q. The connection of this equation with the 
star-triangle relation will be considered in Section 5.1. Note here two important properties 
of the i?- matrix (I3.16p . Introducing 

uj = 1/ql, io'/^ = -1/qo, (3.52) 
where qo is defined in f l2.14p one can show that 

r(A,a;"x) _ f 1 + AV A ^ ^ ~ uj-^^'^X-^xuj^ 



fi{x,z) \l + X-^x^ 
and 



N~l 

Y[f{X,uj^x) = 1 . (3.54) 

j=0 



16 




Figure 3 - The square lattice £ of M rows with L sites per row. Tq is the transfer matrix 
of an odd row, Tq of an even row. Three vertical and two horizontal dotted rapidity lines 
are shown. 

4 Chiral Potts model. 

We define the chiral Potts model in the usual way [6], [15]. Consider the square lattice 
£, drawn diagonally as in Fig. [31 with L sites per row. At each site i there is a spin cr,, 
which takes values 0, . . . , — 1. There is an associated lattice C denoted by dotted lines, 
such that each edge of C passes through a vertex of C 

To each vertical (horizontal) line on the "dotted" lattice C assign a rapidity variable 
p (g). In general they may be different for different lines. In fact a convenient level 
of generality that we shall use here is to allow the vertical rapidities to be alternating 
p,p' ,p,p' , . . ., as indicatecjfl. Each edge of C is assigned a Boltzmann weight depending 
on the two spins adjacent to the edge and on the two rapidities passing through the edge. 
For example, consider a typical SW NE edge of C (with j above i). The spins 
(Ti, aj interact with Boltzmann weight Wpq{ai — aj) (or Wp'q{ai — crj) )■ Similarly, on all 
SE — >• NW edges the spins interact with Boltzmann weight Wpq{ai — aj) (or Wp'q{ai — aj) 
), where again j is above i. The explicit form of the Boltzmann weights will be given 
below. 

Each rapidity variable p (or q) is represented by a four-vector p = {ap, bp, Cp, dp) which 
specify a point on the algebraic curve Ck defined by any two of the following four equations 
(the complimentary pair of equations follows from the other two) 

ttp + k'bp = kdp , k'a^ -\- bp = kc^ , 
kdp ~t~ k Cp — dp , kbp ~t~ k dp — Cp , 

where k"^ + k''^ = 1. The modulus of the curve, k, is considered as a fixed parameter of 
•^Here and below th.e letter q denotes a rapidity variable; it is different from the parameter q in (|2.2p . 
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the model. It is convenient to also use another set of the "p- variables" , 



Xp cip/ dp, Up hpl Cp, Sp dp/ Cp l^'^) 

and 

tp XpUp ttphp/ Cpdp. (4-3) 
In these variables the curve (14. ip reads 

+ = k{l + x^yp, kx^ = l-k's;^, ky^ = l-k's^ (4.4) 
With these definitions the Boltzmann weights have the form 



Wpq{n) -A- COapdq — dpttqU^ -A- UJXp — UJ^Xq 



Wpq{0) cpbq - hpCqUi " Vq - U^Vp 



(4.5) 



P 



Wpqjn) _ -A- dphq - apCqUJ^ _ -A- yg - u^x. 

~ f},bpdq-CpaqU^ -IV^d iYy^_^J^^ 



where n G Z and is a primitive root of unity of degree N . The weights satisfy the 
periodicity conditions Wpq{n + N) = Wpq{n), Wpq{n + N) = Wpq{n). They also satisfy 
the star-triangle relation [6], [26], [16]: 



N-l 



^ Wqr{h - d)Wpr{a - d)W pq{d - c) = RpqrWpq{a - h)W pr{h - c)Wqr{a - c) . (4.6) 



d=0 



for all rapidities p, g, r and all integers (spins) a, 6, c. Here Rpqr is a spin-independent 
function, defined by 

Rpqr fpqfqr/ fpri (4-''') 



where 



and 



N-l 



fv,= \{ {wfqij)IWpq{j)) (4.8) 



j=0 



W^l^{n) = Y,Wpq{a)u:- . (4.9) 

a=0 



It follows from fHTSj) that 



("■) -TT c„6„ - a„d„uj^ ^ ya- UJ^ 



^t?(0) \}^Kcq - dpaquji f}^yp-ujixqSpSq 

(J) 

Note also, that the normalization factor Wpq (0) can also be written in a product form 
by using the identity (2.44) of [45j, namely 

N-l N-l , _ j 

n = WpqiOfN^/' ,..(iv-i)(A^-2)/i2 -Q [t. u^t,} ^4,1) 
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We define row-to-row transfer matrices T and T as in [35]. Let a = (Xi, . . . , ctl be tfie 
spins in tlie lower row of Figure [3l Similarly, let a' = o"^, . . . , be the spins in the next 
row, and a" = a'(, . . . ,a'l those in the row above that. Let T be the A^^ by A^'^ matrix 
with elements 

L 

T.a' = n Wpii^J - (^'jWp'.i^J+i - ^j), (4.12) 
J=l 

similarly, let T be the A^^ by A^'^ matrix with elements 

L 

= n ^pi(^'j - o^p'^i^'j - ^i+i) (4-13) 

J=l 

Let y be a formal variable such that = 1. Define 

N-lN-l N-l 

F{p,q;Y) = A^-^ .;-^iy,,(a)r'' = Y^^iP i-^)"^^ 

a=0 6=0 a=0 ]^^^ 

N~IN~1 N-l 

F{p,q-Y) =N-'Y.Y.^'"''^P'!(^)^' =E^P''(«)^'^ 

a=0 b=0 a=0 

where Wpq (n) is defined similarly to (14. 9 p . It is easy to see that 

F{p, q- u^) = Wp,{n), F{p, q; w^) = Wl[\n) (4.15) 

for any n G Z and that the functions F and F are uniquely determined by the following 
recurrence relations 

F{p,q;Y) _ Spjvg - XpY) F(r, n- 1^ - W (0^ (Aia\ 
Fip,q;u-^Y) - s,iyp-.,Yy ^' " ^-^^^ ^^"^^^ 

Fip,q;Y) _y,-x,SpS,Y p^^^^.^^^^U)^^y 



F{p,q;uj-^Y) Vp-XgSpSgY 

Now choose the following matrix realization of the algebra (I2.32p with q^ = uj~^, 
u;^ = l, 

iX)a,b = 5a,6+l, {Z)a,b = ^"Kh, (4-18) 

where a, 6 G Z^r and 

f 1, a = h (mod A^), 

5a,b = < (4.19) 
I 0, otherwise. 

From and <KW^ it follows that 



N-l 



{F{p,q-Z))a^b = KbWp,{(^), {F{p,q,X))a,i,= Y,Wpg{k)6aMk- (4-20) 



fc=0 



Using these one can equivalently rewrite [16] the star-triangle relation (l4.6p as the following 
relation between F and F 

F{p, q; Z-')F{p, r; X)F{q, r; Z'^) = R-g\F{q, r; X)F{p, r; Z-')F{p, q; X) . (4.21) 
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5 The connection of the discrete quantum Sine-Gordon 
and chiral Potts models 



5.1 The star-triangle relation 

We now want to identify equation (13.491) with the star-triangle relation in the form fl4.2ip . 
Equation fl3.49p contains four continuous parameters A, fi, P, Q, while f l4.2ip involves the 
modulus of the curve (14.ip and three rapidities p, q, r, representing points on that curve. 
We have to establish a relationship between these two sets of parameters. Denote 



(5.i: 



A=l + (e«/A)^, B = l + (e^A)^ 
C = l + (e^//i)^, D = l + (e^/i)^ 

and consider the curve fl4.ip with the modulus 

^ [ACjl -B)- BDjl - C)][BC{1 -A){l-D)~ AB] 

CD[A- B{1- D)][B - A- BC{1- A)] ^ ' ' 

Now choose three rapidities variables p, q and r such that 



UJ"^/"^ jJi~^e^ = XqSqSr/Vr, UJ~^/'^jJie^ = XrSrSq/Uq, 



(5.3) 



where uj and go are related by (I3.52p . One can check that the last four equations are 
consistent with (15. 2p . In fact, the expression (15. 2p is a corollary of (15.30 . This can be 
verified by direct substitution of (15. 3p into (15.10 and (15. 2p with an account of the relations 
(14. 4p . Comparing now ( ]4.10p and the second equation in ( 14.50 with equation ( 13.530 we 
see that the relations (15. 3p allow us to identify the first factors in the left and right sides 
of (I4.2ip with the corresponding factors in ( I3.49P 

F(p, g, Z-^) = r(A, e^^-^), F(g, r, X) = r(/i, e^X), (5.4) 

provided the normalization factors ^^^^(O) and Wqrifi) (which so far were at our disposal) 
are appropriately chosen to satisfy (13.540 . To identify the remaining factors in (I3.49P and 
(14.210 we need four more pairs of relations 

UJ~^/'^{\^)~''e^' = XpSpSr/Vri UJ~^^'^ Xfie^' = Xr-SrSp/Up, 

(^-l/2;\-lgP" — XpSpSq/Uq, LO^^/'^Xe^' = XrSrSp/y, 



(5.5) 



pi 



UJ 



Mysteriously enough, all these relation are corollaries of (15.30 . (I3.5ip and (14.40 . The 
calculations are simple but too long to be presented here. We quote only a few formulae 
important for the derivation of (15.50 . It follows from (15. 3p that 

A' = tq/tp, = tr/tq, (5.6) 

where the t-variables are defined in (14. 3 p and that 

= UXpXq/ypl/q, e^^ = uXqXrslsl/yqi/r. (5.7) 
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Then one can rewrite fl3.5ip as 



e *^ ^Xp^r / VpUri 6 LOXpXrSpS^ I ypUj-i 

e2<9" = uXpXr/ypVr, e^^" = uXpXgslsl/ypyq. 
Similarly to (15.4p we have from (15. 5p 

F{p, r, X) = r(A/i, e^'X), r, Z"!) = r(A/i, e^'Z"!), 

r, Z-i) = f (/i, e«"Z-i), g, X) = r(A, e^"X). 



(5.8) 



(5.9) 



Thus we have proved the equivalence of the twisted Yang-Baxter equation 03.49P with the 
star-triangle relation of the chiral Potts model (I4.2ip . Note that the normalization of the 
weights of the chiral Potts model assumed by the relations (15.40 . (15.90 is such that the 
factor Rpgr, (14. 7p in (I4.2ip is equal to 1, as it follows from (13.540 . 

Let us summarize the results of this section. The four relations (15. 3p allow us to express 
the parameters k, p, q, r entering the star-triangle relation (14. 6 p or (14.210 through A, fi, P 
and Q from the twisted Yang-Baxter equation (13.490 and vice versa. Remarkably, these 
four relations imply (15. 5p enabling us to prove the equivalence of the two Yang-Baxter 
equations (I3.49p and (14. 6p . This equivalence allows us to look at the chiral Potts model 
from a new angle. 

5.2 The evolution operator 

Consider now the evolution operator ( 12.130 . Substituting there the expression for the 
i?-matrix (13.140 we can factorize U into classical and quantum parts [30] 



U = UelUquant (5.10) 

where 

Uci = exp(-7^SV/i) ew{-n^^/h) (1 + 0(/i)), (5.11) 

L-l 

= -Y.H{H.\w2n+k), fc = 0,l, (5.12) 

n=0 

with the function H defined in (I3.47P and (13.150 . Further, 

L-l L-l 

Uquant = W r{K^ , Tcl{w2n)) Y\_ ^('^^ W^2n+l), (5.13) 
n=0 n=0 

rd(a) = fexp(7^iJ)) oexp(H^f ) oa) = limU^J^aUd, (5.14) 

where we used the notation (I2.26p . Let us now remove the constraints (12. 8p and consider 
a special case of the realization (I2.18P when q ^ Qq, q^^ = 1, 

W2n = «2n-^n+l, W2n+1 = a2n+l^r7+l^n+2, n = 0, . . . , L - 1, (5.15) 

such that 

a2n = a, a2n+i = Vn , (5.16) 
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where a and (5 are fixed constants. It follows then from (12.101) that 



r(a2n)=a, r(a2n+i) = Vn, (5.17) 

i.e., the configuration (I5.16P is stationary with respect to the classical evolution. 
Consider the algebraic curve (14.11) with the modulus 



a^/3^) (1 - a^/j- ^^ 
(1 -a^/€2JV)(i _«Nk-2JV) 



and choose two rapidities p and g such that 



uj~^l'^ali<? = XpSpSg/ug, uj-'^/'^aK^ = XgSpSg/yp, 



(5.19) 



Note that the last four equations are consistent with the relation (15.181) and imply the 
latter as their corollary. The expression (15.131) can now be rewritten as 

L L 

Uq^ant = l[F{p,q,X^) F {p, q, Z^^+i) , (5.20) 

n=l n=l 

where we have used (I4.14p . Taking into account (14.201) the matrix elements of this operator 
are 

bi,b2,...,bL ^ ^ 



U quant 



llW^gian-bn) llWpg{bn+i-bn), {5.21] 



J ai,a2,.--,OL 

n=l n=l 



where Wpg{a — h) and Wpg{a — h) are the Boltzmann weights of the chiral Potts model 
defined in (14. 5p . 

Let M > 1 be an integer. Obviously, the trace 

r 1 

^chiral Potts = '^^ Uquant (5.22) 

is the partition function of the chiral Potts model for a (non-diagonal) square lattice of 
size L by M with periodic boundary conditions in both directions (to avoid confusion, 
note that usually the periodic boundary conditions in the chiral Potts model are imposed 
for the diagonal (i.e, 45°-rotated) square lattice, corresponding to the transfer matrix 

6 Concluding remarks 

We have presented a new interpretation of the chiral Potts model where it is arising as 
a particular case of a more general Ising-type model on a square lattice with local spins 
taking N > 2 values at each site. The Boltzmann weights of this model are determined 
by solutions of the classical discrete sine-Gordon model (which is an integrable model of 
classical field theory). In this setting the chiral Potts model corresponds to the simplest 
(constant) solution of the above classical model. It would be interesting to consider more 
general spin models, corresponding to non-trivial solutions of the "background" classical 
field theory. Our construction also sheds some light on the origin of the "non-difference" 
property in the chiral Potts model. Further discussion of this point will be given in [39] . 
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Appendix A. Derivation of Eqs.(13^)-(13^). 



The derivation of the first equation in (13.81) for $ = <l>i goes through the following 
sequence of transformations by using (13.31) and (13. 2p 



$L n = r(A, u)r(A/i, f )r(/i, u) u 

= ur{X, u){c\Xfi + v){q + A/if )^V(A/i, v)r{fi, u) 



, . 1 + qXu\ f , qA + f 
u qA/i + V— M + -^A^"" 



X + qu 



q + Af 



r(A, u)r{Xfi, v)r{fi, u) 



ug{u,v,q)^L 



where the function g is the same as in (13.41) . 

For $ = the sequence of transformation is longer 



(A.l) 



u = r{fi, v)r{Xfi, f )r(A, f ) u 

= r(/i, f ) ur{X^, m) '^^"^ r(A, f ) 

= r(/i, f ) u (qA(A/x + qu) + v{l + qXfiu))x 

(q(A/i + qu) + Af (1 + qXfiu))^^{r{Xfi, u)r{X, v] 

r(/i, f )(qA(A/i + qu) + q^f (1 + qXfiu)) u 

(qAV + <fv + qA M(q/i + v))r{fi, v) u 

q/i + f 



u (qA^/i + q^A/iM + v + qXvu 



-j>r(/i,' 



q + /if 
(q^ + AVq-M) '^R 



. . . M A(q/x + v)+u 

J V q + /if 

r ^ | q + /if _^ _^ _^ q^A^ff )"^$ 
I J q/i + f 

f5((f,f, q)$R. 



(A.2) 



The second equation in 03.80 is obtained from (lA.ip and (1A.2P with the help of the 
automorphims ( 13.60 . 
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